Anisotropic Exchange Interactions of Spin-Orbit-Integrated States in Sr2lr04 
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We present a microscopic model for the anisotropic exchange interactions in Sr2lr04. A direct construction 
of Wannier functions from first-principles calculations proves the jeff=l/2 character of the spin-orbit integrated 
states at the Fermi level. An effective jeff-spin Hamiltonian explains the observed weak ferromagnetism and 
anisotropy of antiferromagnetically ordered magnetic state, which arise naturally from the jeff=l/2 state with a 
rotation of IrOe octahedra. It is suggested that Sr2lr04 is a unique class of materials with effective exchange 
interactions in the spin-orbital Hilbert space. 



I. INTRODUCTION 

Many of transition metal oxides (TMOs) are antiferromag- 
netic insulators. The simplest model for such Mott insula- 
tors is the Hubbard model Hamiltonian/ which gives rise 
to an effective exchange term called superexchange interac- 
tion at half-filling. When orbital degrees of freedom are in- 
volved, a variety of exchange interactions can occur for a 
given ionic configuration with different crystal structures. In 
the case of colossal magneto-resistance manganese oxides, for 
instance, the superexchange interaction with orbital degener- 
acy determines complex spin and orbital orderings and, when 
doped, degenerate Cg orbitals coupled to the lattice via Jahn- 
Teller interactions become an essential part of the double ex- 
change physics^. Sometimes the orbital degrees of freedom 
via spin-orbit (SO) coupling are responsible for the magnetic 
anisotropy bound to the crystal environment. When there ex- 
ists an orbital degeneracy, SO coupling may become a dom- 
inant term so that the effective Hamiltonian should involve 
the full spin-orbital Hilbert space where the ground state must 
comply with the intersite spin and orbital correlations^ . A pos- 
sible dynamic interference between the spin and orbital space 
was suggested in vanadates^. There was a report of a large 
spin-orbital fluctuations in Mott insulators with t2g orbital de- 
generacy as a manifestation of quantum entanglement of spin 
and orbital variable^. 

Recently we have shown that the electron correlation ef- 
fect combined with strong spin-orbit (SO) interactions is re- 
sponsible for the observed insulating behavior oi 5d TMO 
Sr2lr04^. While SO coupling has been considered as a mi- 
nor perturbation in the description of magnetism^, the amount 
of SO interactions m 5d elements including Ir, for example, 
is an order of magnitude larger than in the ?>d TMO system^. 
Thus the SO coupling is expected to play a significant role in 
the electronic and magnetic properties of 5d TMO systems. 
Indeed the manifestation of a novel jeff=l/2 Mott ground 
state in Sr2lr04 was revealed by angle resolved photoemis- 
sion spectroscopy, optical conductivity, x-ray absorption spec- 
troscopy measurements, and first-principles electronic struc- 
ture calculations^. Further investigations of the electronic 
structures of the Sr^+iIrnOsn+i (n = 1, 2, and oo) series 



demonstrated a Mott insulator-metal transition with a change 
of bandwidth as n increased. The ground state of 5d TMO 
Sr2lr04 is a Mott insulator in the strong spin-orbit coupling 
limit. In addition, Sr2lr04 exhibits unus ual weak fer romag- 
netism with reduced Ir magnetic moment j i Q | 1 1 1 1 2 | 1 3 1 i H 
derstand such unusual magnetic properties of Sr2lr04, it is 
necessary to take account of the spin-orbit integrated jeff=l/2 
state. 

In this paper, we introduce a prototype model of spin- 
orbit-integrated magnetism realized in Sr2lr04. From a tight- 
binding analysis based on first-principles calculations, we 
show that the jeff=l/2 character of the spin-orbit-integrated 
state remains robust even in the presence of on-site Coulomb 
interactions. A direct construction of Wannier functions from 
first-principles calculations proves the jeff=l/2 character at 
the Fermi level. An effective exchange Hamiltonian with not 
5'=l/2 but jeff=l/2 is obtained starting from a jeff=l/2 Hub- 
bard model. The origin of anisotropic magnetic exchange in- 
teractions are discussed in connection with an extraordinary 
character of the ground state. The presence of spin-orbit- 
integrated state with strong SO interactions in Sr2lr04 can 
make 5d Ir-oxides a unique class of materials for the study of 
effective exchange interactions in the full spin-orbital Hilbert 
space. 



II. SPIN-ORIBT-INTEGRATED ELECTRONIC STATES 

A. LDA+SO+[/ Band Structure 

Since both on-site Coulomb interactions (U) and SO cou- 
plings are expected to be important in the description of Ir 
5d states, we examined the effect of on-site U and SO cou- 
plings separately and simultaneously on the electronic struc- 
ture of Sr2lr04. To identify the role of each term and the 
interplay between them, we carried out density-functional- 
theory (DFT) calculations within the local-density approxi- 
mation (LDA), LDA including the SO coupling (LDA+SO), 
and LDA+/7 including the SO coupling (LDA+S0+/7), re- 
spectively. We calculated total energies and electronic band 
structures of Sr2lr04 for the structural parameters as obtained 
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FIG. 1: (Color online) Calculated Wannier functions of the |jeff = 
1/2, +1/2) state: (a) the real part of up-spin | j) component, (b) the 
real part of the down- spin | |) component, and (c) the imaginary part 
of the down- spin | I) component. Blue (dark gray) and yellow (light 
gray) colors in the Wannier function represent negative and positive 
values respectively. The inset shows the LDA-i-SO-i-[/ band structure 
near ^f=0 eV, emphasizing the "upper" and "lower" Hubbard band 
of the I jeff = 1/2, m) band above and below respectively. 



from the neutron powder diffraction data at 10 K^^, which 
has a K2NiF4-type layered perovskite structure with the sym- 
metry of the space group Mi/acd reduced from M/mmm, 
where IrOe octahedra are rotated by about 11° around the 
c-axis of the unit cell. For the calculations, we used the 
DFT code, OpenM^^PI, based on the linear-combination-of- 
pseudo-atomic-orbitals method^ where both the LDA+[/ 
method^^ and the SO couplings were included via a relativis- 
tic i -dependent pseudo -potential scheme in the non-collinear 
DFT formalisnPEOHa. Double valence and single polariza- 
tion orbitals were used as a basis set, which were generated 
by a confinement potential scheme with cutoff radii of 8.0, 
7.0 and 5.0 a.u. for Sr, Ir, and O atoms respectively. We used 
a (6x6x4) k-point grid for the k-space integration. 

Calculated LDA, LDA+SO, and LDA+S0+/7 band struc- 
tures were presented in the previous work^ where the re- 
sults of the LDA+S0+/7 band structures of Sr2lr04 are well 
compared with those of angle resolved photoemission spec- 
troscopy experiments. The LDA bands near the Fermi level 
(Ef), as shown in Fig. 2(a) of Ref. 1 8 1, are almost identical to 
those of Sr2Rh04,^^ which can be expected from the same 
configuration of Rh^+ and Ir'^^ and the same structural dis- 
tortions, i.e., the rotations of RhOe and IrOe octahedra. The 
hybridization of dxy and d^'^-y^ due to the rotation of IrOe oc- 
tahedra pushes the d^y band below E^, similarly to the case 
of Sr2Rh04. Indeed the LDA Fermi surface of Sr2lr04 was 
found to be basically the same as that of Sr2Rh04^^. 

In the LDA band structure, the contribution of dxy com- 
ponents above E^ are strongly suppressed relative to those 
of dyz and dzx states, whereas the Ir 5d bands ranging from 
— 2.5 eV to 0.5 eV are still dominated by the t2g orbitals with a 



TABLE 1: Coefficients of the Wannier functions illustrated in Fig.[T] 
Only the coefficients from the center Ir are listed. 
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small admixture ofdx'2-y'2. On the other hand, however, when 
the SO coupling is included, a significant change of the wave 
function character occurs so that all three t2g orbital com- 
ponents are almost equally distributed in the LDA+SO band 
structure. This change arises from the SO interactions acting 
on the t2g manifold, which mixes up the dxy, dyz, and dzx or- 
bitals. This qualitative change of the wave function character 
is an essence of the SO coupling action, which is related to the 
novel nature of the SO-integrated insulating ground state. 

As shown in the inset of Fig. [T] an effective U = 2 oV 
opens up a gap in the LDA+S0+/7 band structure and gives 
rise to the non-dispersive and parallel features of "upper" and 
"lower" Hubbard bands of the SO-integrated states, which are 
in excellent agreement with experimental observations^. It is 
remarkable to obtain an insulating ground state for the inter- 
mediate value of /7, which is smaller than the band width of 
the t2g manifold and the conventional U values of 3d TMOs. 
On the other hand, however, when we performed LDA+/7 
calculations without the SO coupling, the on-site Coulomb 
interaction became ineffective due to the three-fold degener- 
acy of the Ir t2g manifold crossing Ejf ^^. Unless the degen- 
eracy is broken, each band remains partially filled being far 
from the Mott instability. In our non-collinear DFT calcu- 
lations, the SO coupling terms were solved in a completely 
non-perturbative way, whereas the Coulomb correlation effect 
were treated via the LDA+/7 method. 

In order to examine the nature of the SO-integrated state, 
we constructed Wannier functions which can identify the or- 
bital shape and bonding character of the "upper" and "lower" 
Hubbard bands of the SO-integrated states, as shown in the 
inset of Fig. [T] The Wannier functions were calculated for 
the t2g manifold by employing the projection scheme^"^. The 
Wannier function illustrated in Fig. [T] corresponds to the sin- 
gle band above Ep of LDA+SO+[/ bands. As Hsted in Tabled) 
the overall shape of the calculated Wannier function of Fig. [l] 
matches closely to the ideal jeff=l/2 state: 

beff = l^^l) = ^^ [\dxy)\±) ± {\dy,) ± i\d,.)) It)] , 

(1) 

where |±) represent for the up-spin | |) and down- spin | |) 
states respectively. The agreement of its orbital components 
and their relative phases between the ideal state and the calcu- 
lated Wannier function is another proof of the SO-integrated 
j^^=l/2 state. Here, for the sake of simplicity in the presenta- 
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tion, we chose a self-consistent solution with the spin quanti- 
zation axis parallel to the z-axis. 



B. Tight-binding model 

The physics of the LDA+SO+[/ results can be captured by 
a multi-band Hubbard model for the t2g bands including the 
SO coupling term. The tight-binding (TB) bands for the t2g 
manifold can be described by 
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(2) 

where (ij) runs over the nearest neighbor pairs of sites i 
and j in the two-dimensional square lattice consisting of Ir 
ions, a and (3 are indices for t2g orbitals, i.e., {d^yj dyz^dzx}, 
t^^ a hopping integral between \ia) and At a tetrago- 
nal crystal field splitting, i.e., an on-site energy difference of 
the dxy orbital relative to dyz and dzx, and Aso the SO cou- 
pling parameter. In a simple square lattice of Ir ions, be- 
comes a non-zero constant to only for (a, P) = {dxy, dxy) = 
{dzx-, dzx) with j = z + X and so on. 

Starting from a set of {{uiraa)} as mean-field parame- 
ters, we could obtain a self-consistent mean-field Hamilto- 
nian within the t2g subspace by Ht2g = Xlk^k^(^)^k 
where Ck has 12 components of {ckraak = A^B;a = 
dxy J dyz, dzx] o- =T, i}. Here the site indices r = A^B are 
for the two inequivalent Ir sites. By choosing the basis in or- 
der of {cAdxy^,CAdy,i,CAd,x^, [A B], [t^j]), we can find 
a block-diagonal 12x 12 T(k) matrix: 
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Here the hopping integrals contribute to 0(k): 
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where the non-zero hopping terms of to 
td^^ lead to the dispersions 71k 

for d. 



cos ^ cos 

cos ^^1^, and 73k = cos ^^2^ x^x ^^y, u.yz, 
respectively. The on-site Coulomb interaction U and the SO 
coupling Aso contribute to the diagonal term: 
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(5) 



and D[j is a time-reversal partner of DJ" . 

Despite of a large cubic crystal field splitting due to Ac « 5 
there is a significant hybridization of 
y2 due to the rotation of IrOg octahedra. In order 
to describe both LDA and LDA+SO band structures properly, 
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FIG. 2: (Color online) (a) Tight-binding band structure with (Aso, 
U)={QA eV, 2.0 eV), which is well compared with the LDA+SO+[/ 
band structure of Fig. [ijc), (b) the decomposition of the "upper" 
Hubbard band wavefunction (marked by a thin (red) solid line in (a)) 
into the {| jeff, ^j)} basis, where the dotted line represents a ||, |) 
component with [7=0, and (c) a schematic energy diagram of the t2g 
manifold in the atomic limit. Due to the large crystal field splitting 
by Ac, the t2g levels can be mapped into the effective /eff = 1 states 
where the tetragonal field splitting by At is relatively insignificant. 



the contribution of the dx2-y2 admixture is necessary to be 
included as a k-dependent energy Ae\^dxy the dxy band: 
Ae\idxy — ^i7ik- ^^^t fit to the LDA bands was obtained 
by a set of parameters: At= 0.15eV, to=035eV, and ei = 
-1.5eV. 

The solutions of our TB model including both SO coupling 
Aso and on-site Coulomb interaction U with different sets of 
parameters (Aso, ^)=(0, 0), (0.4 eV, 0) and (0.4 eV, 2.0 eV) 
reproduce well the t2g manifold of the LDA, LDA+SO, and 
LDA+S0+/7 bands, respectively. The self-consistent solution 
for (Aso, U)=(OA eV, 2.0 eV) is shown in Fig.[2];a), which 
corresponds to the LDA+S0+/7 bands of the inset of Fig. [T] 

In addition to the large crystal field splitting between Cg 
and t2g, the t2g manifold splits further into doubly degenerate 
jgff=l/2 and quadruply degenerate jeff=3/2 states due to SO 
coupling. The small tetragonal crystal field does not affect 
this configuration. A schematic energy level diagram shown 
in Fig.[2jc) was confirmed by the LDA and LDA+SO energy 
levels at the X point, ^ where the off-diagonal hopping ma- 
trix 0(k) in the TB model becomes zero. Even though the 
non-zero hopping terms away from X point may disturb the 
atomic picture, the SO coupling retains the anti-crossing be- 
tween those levels which transform according to the same ir- 
reducible representatioiP^. Consequently, the effective band- 
width of the half-filled jeff=l/2 band becomes smaller than 
the modest value of on-site U. The decomposition of the "up- 
per" Hubbard band wavefunction into the {Ijeff^^j)} basis 
clearly demonstrates the robustness of its jeff=l/2 character 
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as shown in Fig.|2jb), whereas the jeff=l/2 weight for = 2 
eV is sHghtly reduced from that ofU = 0. Therefore it is 
reasonable to consider an effective Hamiltonian based on the 
j^^=l/2 single-band Hubbard model instead of the conven- 
tional S=l/2 model: 

^= "^ini'Arndjni' ^U^ndi+l/2ndi-l/2. (6) 

{ij)mm' i 

where dim represents for the |jeff = 1/2, m) state at the site 

i with m, m' ±1/2 and Udim Am^irn- ^mm' ^"^^ ^ 

effective hopping and on-site interaction parameters respec- 
tively. 

III. ANISOTROPIC EXCHANGE INTERACTIONS 

A. Effective Exchange Hamiltonian 

The Jeff =1/2 single-band Hubbard model has an interesting 
feature in t^^/ , which originates from a peculiar nature of the 
spin-orbit integrated state in Sr2lr04. In the strong SO cou- 
pling limit, the orbital wavef unctions of the |jeff = 1/2, m) 
state of Eq. ([T]) consists of the cubic harmonics with respect 
to the local coordinate axes. The rotation of the IrOe oc- 
tahedron results in a rotation of the |jeff = 1/2, m) state 
at each site i, thereby generating a spin-dependent hopping 
term. In Sr2lr04, where the IrOe octahedron is rotated by a 
angle ^ 11° about the c-axis, the effective hopping ma- 
trix t^^/ can be represented in terms of Pauli matrices by 
= fol + ^ficr^ where fo and ti for {ij) = x or y become 

tt/^ = ^cos^(2cos^^-l) (7) 

ti^^ = ^sme{2sm^e-l). (8) 
o 

At half-filling, an effective jeff-spin Hamiltonian can be de- 
rived from the jeff=l/2 single-band Hubbard model of Eq. ([6]): 

^spin = ^ [lo^i • Jj + hJziJzj + ^ij ' Ji X Jj] (9) 

{ij) 

where Iq = 4@ -P^)/U,h = Stf/U, and D^^- = D^z with 
Dz = Stoti/U. The first term is a conventional Heisenberg 
form of superexchange with the coupling constant /q. The 
second and third terms are pseudo-dipolar and Dzyaloshinkii- 
Moriya (DM) antisymmetric exchange interactions, which 
originate from the pure imaginary hopping matrix element iti 
between the neighboring | jeff = 1/2, m) states of rotated IrOe 
octahedra. 



B. Comparison with LDA+SO+[/ Results 

From our LDA+S0+/7 calculations, the magnetic configu- 
ration of the insulating ground state was determined to be a 




FIG. 3: (Color online) (a) Magnetic configuration and (b) DM vec- 
tors of the calculated LDA+SO+[/ ground state of Sr2lr04. Blue 
arrows in (a) represent a non-collinear ordering of the local Ir mo- 
ments, consisting of both spin and orbital components, in a canted 
AFM configuration. The DM vectors. Dab and Dbc, in (b) are 
aligned along the c-axis with alternating signs and consistent with 
the Dzyaloshinkii-Moriya (DM) rule. 



canted antiferromagnetic (AFM) state with the a6-plane as an 
easy plane. We found no preferred direction within the ab- 
plane. As illustrated in Fig.[3ja), there are two inequivalent Ir 
sites, i.e., Irl and Ir2 within the a/2 x a/2 unit cell. It is found 
that the magnetic moment at each Ir site is 0.36 /iB and both 
spin (0.10 /ie) and orbital (0.26 /ie) moments are parallel to 
each other. In addition, AFM moments are canted with the 
canting moment 0.063 /iB, which is comparable to the single 
crystal measurement^^. 

According to the rule by Dzyaloshinkii and Moriya^, the 
direction of the vector D^^ in Sr2lr04 should point to the c- 
axis due to a mirror plane containing Irl-0-Ir2, as illustrated 
in Fig. |3|b). The directions of T>ij can be represented by 
Dab = — Dbc = (0, 0, dc), when considering the inversion 
symmetry at the site B, which gives the consistent results as 
the jeff=l/2 Hamiltonian of Eq. From the LDA+S0+/7 
calculations, it is concluded that the DM interaction is respon- 
sible for the magnetic anisotropy of Sr2lr04 with the a6-plane 
as an easy plane but isotropic within the a6-plane, whereas the 
single-ion anisotropy term has a negligible contribution. Con- 
trary to the La2Cu04,^^ which has no single-ion anisotropy 
due to the 6'=l/2 ground state, the absence of the a6-plane 
anisotropy in Sr2lr04 is attributed to the tetragonal symmetry. 

From the effective exchange Hamiltonian of Eq. (|9]), the 
ratio of Dz/Io, which determines the spin canting angle, be- 
comes \Dz/Io\ ~ tan 2^ for small 6. In the strong SO cou- 
pling limit, the canting angle increases close to the rotation 
angle of IrOe octahedra. We can estimate the magnitude of 
Dij to be |D| ^ 3.8 meV assuming the intersite superex- 
change interaction J « 10 meV. This enormous DM interac- 
tion may well be related to the peculiar nature of the jeff=l/2 
state. Contrary to the 5*= 1/2 counterpart of La2Cu04,^^ the 
j^^=l/2 state has an open- shell of the / = 1 orbital where the 
non-perturbative ground state of jeff=l/2 spin-orbit coupled 
state contribute to the DM term. Although the small magnetic 
moment of Ir observed in experiments was attributed to the 
effective moment the jeff=l/2 state, one can expect possible 
contributions from the j=l/2 quantum- fluctuation. Neverthe- 
less, since jeff=l/2 state is an eigenstate of the fictitious angu- 
lar moment Jeff = Leff + S = —L + S, the orbital contribution 
to the magnetic moment needs a careful interpretation^. 
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IV. CONCLUSIONS 

In summary, we presented the effective jeff-spin model 
Hamiltonian for Sr2lr04. The strong SO interaction combined 
with the large crystal field splitting in 5d TMOs introduces a 
unique form of the spin-orbit integrated band state at E^, lead- 
ing to an effective insulating ground state of jeff=l/2 quantum 
magnet. The observed weak ferromagnetism is understood by 
the DM anisotropic exchange interaction where the effective 
exchange interactions arise from the full spin-orbital Hilbert 
space. We hope that our prototype model of the spin-orbit 
integrated magnetism is useful for the study of various spin- 
orbit entangled physics. By taking an analogy of the high Tc 



superconductors as a doped S=l/2 quantum magnet, it will 
be interesting to observe a doped j=l/2 quantum magnet as a 
spin-orbit integrated correlated electron system. 
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